Typeset with jpsj2.cls <ver.l.2> Full Paper 

Evaluation of Dynamic spin structure factor 
for the spin-1/2 XXZ chain in a magnetic field 
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Transition rates and dynamic spin structure factor at zero temperature for the spin-1/2 
XXZ chain at critical regime in a magnetic field are numerically evaluated in terms of the 
exact determinant representations for the form factors and norms of the Bethe eigenstates. 
We have seen that the transition rates converges toward the constant function with the value 
1 in the limit A — * 0. The observed critical exponent of the singularity at the lower boundary 
is compared with the one predicted from the comformal field theory. We confirm that they 
are in good agreement. Further we have discovered that a small peak emerges near the upper 
boundary in the line shape of S(q,oj) for < A < 1. 
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1. Introduction 

Since Bethe 1 in 1931 formulated the method (Bethe ansatz) for obtaining the exact eigen- 
values and eigenvectors of the one-dimensional spin-1/2 Heisenberg model, it has been ex- 
tensively generalized to study the excitations and the thermodynamics of the model exactly. 2 
However, from the point of view of the dynamics of Heisenberg antiferromagnetic chains, only 
approximate results had been obtained until recently, partly due to the lack of the exact 
expression for the form factors. 

Heisenberg antiferromagnetic chains are realized in quasi-one-dimensional magnetic insu- 
lators, such as, Cu(C4H4N2)(N03)2, 3 and KCUF3. 4 Dynamic spin structure factor S(q,oj), 
which is the Fourier transform of dynamical correlation functions, is of considerable impor- 
tance, since they are directly comparable with inelastic neutron scattering experiments of 
these quasi one-dimensional substances. 6 

In 1967 Niemeijer obtained an exact expression of S(q, u) for the spin-1/2 XY model at any 
temperature. 5 It is the special case (A = 0) of XXZ model, where all the correlation functions 
can be calculated by Jordan- Wigner transformation. Unfortunately this exact calculation can 
not be extended for general case of A 7^ 0, and only approximate evaluations of S(q,ui) had 
been attempted. 7 It was worth noting, therefore, that an analytical result for the two-spinon 
dynamic spin structure factor for the spin-1/2 massive XXZ model (A > 1) in zero magnetic 
field was obtained 8-10 from the multiple integral representation for form factors, 11 which is 
based on the infinite dimensional symmetries of the quantum affine algebra U q (sl{2, C)). 
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However, this analytical calculation is not applicable to the case with nonzero magnetic 
field. Karbach et al. classified the dynamically dominant excitations of the spin-1/2 XXX 
Heisenberg model (A = 1) in a magnetic field. They subsequently calculated their transition 
rates pertaining to dynamic spin structure factor, directly from the Bethe wave functions for 
finite chains with up to N = 32 sites. 12 ' 13 Although solving the Bethe ansatz equations in 
its own is possible for N ~ 10 3 , constructing the Bethe wave functions is limited to N ~ 32. 
The problem of this direct calculation lies in the evaluation of the sum over the r! magnon 
permutations in the coefficients of the coordinate Bethe wave functions, where r is the number 
of the down spins. 

Quite recently, this problem was solved by the novel work of Kitanine et al., who de- 
rived the exact determinant representation for the form factors of local spin operators for 
arbitrary finite systems. 1 In the framework of the algebraic Bethe ansatz, 15 they have found 
that any local operators can be written as the elements of the quantum monodromy matrix. 
Consequently the calculation of form factors reduces to that of scalar products between a 
Bethe state and an arbitrary state, which has a determinant representation. Then it enables 
us to calculate transition rates directly from the solution of Bethe ansatz equations, without 
constructing the Bethe wave functions. 

Biegel et al. calculated dynamic spin structure factors for chains with the order N ~ 10 3 
using this determinant representation in the case of XXX Heisenberg chain in a magnetic 
field 16 and in the case of XXZ Heisenberg chain at critical regime in zero magnetic field. 17,18 
As much as these calculations are still for the finite systems, we can deal with much larger 
system sizes than in the past. 

In this paper, as an extension of these works, we shall generalize the work by Biegel et al. 
and evaluate dynamic spin structure factors for the spin-1/2 XXZ chain in a magnetic field. 
Under the existance of a magnetic field, only the XXX case has been dealt with so far. 

2. Bethe ansatz 

Let us consider the one-dimensional spin-1/2 XXZ model with periodic boundary condi- 
tions in a magnetic field 

N r i i N 

H = J^2\ S x n S x n+l + SlSl +1 + A(S*S z n+1 - -) - h Y, S z n , (1) 

n=l J n=l 

where J > is the coupling constant, S x ' y ' z = —a x ' y ' z , A is the anisotropy parameter and h 
is the external magnetic field applied to the positive direction of z-axis. a x,y,z represents the 
standard Pauli matrices. 

It is exactly solved by Bethe ansatz method. Bethe eigenstates with r down spins are 
constructed by a set of rapidities {zi, . . . , z r }, which is a solution of Bethe ansatz equations 

A^tan _1 [cot ^tanhzj] = nli + YZ , { tan -1 [cot 7 tanh(zj — Zj)], 
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i = l,...,r, (2) 

where the parameter 7 is related to the anisotropy A as 7 = cos -1 A. Ii are Bethe quantum 
numbers, which have integer values for odd r and half odd integer values for even r. The total 
momentum and energy eigenvalues are given by 

2tt r 

k = 7TT ~ ~N ^ Ji ' ^ 

i=i 

E = J± " Sm27 h(^-r). (3b) 

^ cosh 2zi - cos 7 V 2 / 

Every solution of Bethe ansatz equations is uniquely determined by Bethe quatum numbers 
I\ < ■ ■ ■ < I r , which provide us the classification of the excitations. The ground state \G) at 
magnetization < M z < N/2 is specified by the set of r = N/2 — M z Bethe quantum numbers 
h = -N/4 + M z /2 + i - 1/2. 

3. Dynamically dominant excitations 

In neutron scattering experiment at sufficiently low temperature, the cross section is con- 
sidered to be proportional to the dynamic spin structure factor S(q, uj) at zero temperature, 
which is defined by the space-time Fourier transform of the dynamical correlation function: 

S^u) = / dte**{S£(t)Si,{0)) 



n,n' 
00 



00 



dte iujt {S%(t)Srt{$)) 
= 2vr^|(G|^|A)| 2 ^-^ A ), (4) 

A 

where (n,p) = (z, z), (+,—),(—, +), lu\ = E\ — Eq. Eq and E\ are the energy eigenvalues 
of the ground state \G) and one of the excited states |A), respectively. The spin fluctuation 
operator Sq is defined by 

1 N 

s $ = -mT, eiqns £> v = z,+,-, (5) 



N 1 

n=l 

for wave numbers q = 2ttI/N, I = 1, . . . , N. 

In the following, we will focus on the case for the parallel spin fluctuations (fi, Jx) = (z, z), 
q = it/2, and for half the saturation magnetic field M z = N/4, r = N/2 - M z = N/4. In 
the case of XXX model (A = 1), it was established that the spectral weight of S zz {q, uj) is 
dominated by the set of collective excitations, called " psinon(^f)-antipsinon('^*) v excitations, 
which is denoted by tyfy*. 12 Their Bethe quantum numbers are shown in figure 3. The Ii are 
given by the position of small black circles in each row. The positions of large circles represent 
the Ii vacancies. The first row represents the ground state \G), whose Bethe quantum numbers 
are uniformly cofigured at the center. 
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Fig. 1. Configurations of Bethe quantum numbers of psinon-antipsinon excitations for N = 16, M z = 
A,r — 4. Ii arc represented by the positions of small black circles. Large circles represent 7j 
vacancies. The first row represents the ground state \G). 



When one of the Bethe quantum numbers of the ground state is moved to the sea of 
Ii vacancies, it becomes antipsinon. Then it makes a hole in the Ii configurations, which 
represents the psinon. In this way, we have a set of excitations with two parameters, which is 
called psinon-antipsinon excitations. We label the >I"I'* states with the integer parameter m. 
The state \m) has the Bethe quantum number I\ = —r/2 — m + 1/2. 

The relative contribution of psinon-antipsinon excitations is determined by the ratio of 
the integrated intensity 

/OO J . 
^ -S zz (g,u) = £ |(G|^|A)| 2 = (G\S*Si q \G) (6) 
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Tabic I. Relative contributions of psinon-antipsinon excitations for N = 24, M z = r = 6, q = n/2. 



A 


0/2) 


S zz (ir/2) Rek 


itive contribution(%) 


-0.4 


0.228545 


0.232882 


98.1378 


-0.3 


0.234930 


0.237182 


99.0504 


-0.1 


0.245605 


0.245782 


99.9282 


-0.01 


0.249581 


0.249582 


99.9994 


0.01 


0.250415 


0.250417 


99.9994 


0.1 


0.253984 


0.254114 


99.9489 


0.3 


0.260983 


0.261918 


99.6431 


0.4 


0.264076 


0.265574 


99.4360 


0.5 


0.266954 


0.269056 


99.2189 


0.7 


0.272230 


0.275491 


98.8161 


0.9 


0.277048 


0.281247 


98.5068 


1 


0.279318 


0.283887 


98.3907 


1.1 


0.281503 


0.286378 


98.2977 


1.5 


0.289411 


0.295027 


98.0965 


2.0 


0.297553 


0.303452 


98.0559 



and fyty* contribution 

= £ \(G\S*\W*)\ 2 = ^\(G\S z q \m)\ 2 . (7) 

Relative contributions for various values of A are shown in table I. We see they are more than 
98% for general A. Although the relative contribution actually decreases as the system size iV 
increases, it was shown in Ref.[ll] that the relative contribution for the XXX model (A = 1) 
is more than 93% from an extrapolation of the data for N = 12, 16, 20, 24, 28, 32. Also we find 
that the relative contribution is monotonously increasing when we bring down the value of A 
from 1 to 0. Hence we can conclude that the fyfy* excitation is dynamically dominant in the 
whole region < A < 1. 

4. Transition rates and Dynamic spin structure factor 

We numerically evaluate the dynamic spin structure factor from the contribution of ^E'^'* 
excitations: 

S?** (9, = 2vr £ | (G\ S z q | m) \ 2 5(co - co rn ) . (8) 

m 

In the thermodynamic limit N — > oo, it can be represented by the product of the transition 

2ir 1 

rates M zz (q,u m ) = N\(G\S q \m}\ 2 and the density of states D zz (q,to m ) = — . 13 
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From the determinant representations for the form factors and norms of Bethe eigen- 
states, 14 we have the following formula for transition rates: 



M zz (q,uj) =— — 



J] I smh(z° - % | 2 



II I sinh(^ - — ) 



n 

j<k 



sinh 2 ^ 1 ? 



z%) + sin 2 7 | 



n 

\a</3 



sinh (z a — zp) + sin 7 



det(if - 2P) 



idet #({*?})! |det<2>({z,})| 



(9) 



where the matrix elements of H, P and are given by 

sin 27 



sinh 2 (2; a - z b ) + sin 2 7 
sin 7 



N 



sin 27 



sinh z a + sin ^ sinh (z a — z^) + sin 7 



(a = 6) 



P({z?}, {**})«* 



]^[ sinh(z° - Zf, - ry) - 



sinh(z° - z 6 ) 

, n ~ JJsinh(z fe -z b -ij). 

smh z a + sm - fc=1 



sinh^ 



sinh(z fe + ^) 



A? 



| J sinh(z° - z b + 27) 



(10) 



{^J, . . . , z^} and {zi,...,z r } are the rapidities corresponding to the ground state |G) and 
one of the tyfy* states |m), respectively. To obtain the desired transition rates, we first solve 
the Bethe ansatz equations (2) for given Bethe quantum numbers Ii, and then substitute the 
solution to the formula (9). 

We rewrite the Bethe quantum numbers Ij of the Hfty* state \m) at q = ir/2: 

= - r /2 - m + 1/2 
Ii = -r/2 + i - 3/2 (2<i<r-m + l) 

(11) 

Ii = -r/2 + i- 1/2 (r-m + 2<i<r) 
for 1 < m < r 

When m increases and approaches toward r, the first Bethe quantum number I\ gets away 
from the other Ii configurations, and the transition rate monotonously decreases. At the same 
time it is getting more difficult to find the solution of the first rapidity z\ . 

In figure 2(b), we plot transition rates M zz (it/2,uj) for N = 4096. Also shown in the inset 
(a) are the excitation spectrum of states for each A. Note that in the XXX case (A = 1), 
the transition rates M zz (it/2,uj) can be observed as a smooth function of lo, which has a 
singularity at the lower boundary, and converges toward zero at the upper boundary. 16 
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Fig. 2. (a)Encrgy versus momentum of the psinon-antipsinon excitations at M z — N/4 for 7V=256. 
(b)Transition rates between the ground state and the psinon-antipsinon excitations at q = n/2, 
M z = N/4 for Ar=4096. 
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Table II. Critical Exponent of the singularity at the lower boundary 



A 


1 


0.9 


0.7 


0.5 


n z 


4.38286228 


4.39050106 


4.40585570 


4.42060672 


2irv 


5.72474288 


5.64334684 


5.45713287 


5.23326627 


2d z 


1.53119969 


1.55599192 


1.61471447 


1.68942549 


Vz 


1.5313 


1.5565 


1.6163 


1.6923 



A 


0.4 


0.3 


0.2 


0.1 


0.01 


Q z 


4.42734708 


4.43333735 


4.43823777 


4.44160837 


4.44286904 


2nv 


5.10437806 


4.96243131 


4.80595322 


4.63334574 


4.46278471 


26 z 


1.73472538 


1.78676018 


1.84697501 


1.91723589 


1.99107478 


Vz 


1.7384 


1.7910 


1.8516 


1.9213 


1.9918 



In the XXZ case at the critical regime < A < 1, it can be observed that when we bring 
close to A = starting from A = 1, the singularity at the lower boundary has been weakened 
and at last the function M zz {ir/2,u) converges to the constant function with the value 1. 

The critical exponent of the singularity at the lower boundary can be compared with 
the prediction of conformal field theory. The observed critical exponent i] z in table II is a fit 
M zz (tt/2, oj) = a+uf 1 *" 2 from the data of the lowest and the second lowest excitations m = 1, 2. 
Conformal field theory predicts that it is equal to the scaled energy gap 26 z = Q z /(irv). 19 Q z 
and the spin- wave velocity v are defined by the scaled lowest excitation energy at q = tt/2 
and q = 2tt/N in the thermodynamic limit: 

Q z = lim Noj(m = 1, q = tt/2) 

(12) 

2ttv = lim Noj{m = l,q = 2ir/N). 

N^oo 

They are shown for < A < 1 in table II. We can see that the observed singularity and the 
prediction from conformal field theory agree quite well. 

Next we consider the dynamic spin structure factor S zz (tt/2, oj), which is obtained by the 
product of transition rates M zz (ir/2, oj) and density of states D zz (tt/2, oj). Density of states of 
tyty* excitations are shown in figure 3(a). It begins with flat function from the lower boundary, 
and gradually increases toward the singularity at the upper boundary. 

The spectral-weight distribution S zz (tt/2,lo) is shown in figure 3(b). In the XXX case 
A = 1, it has double singularities at the lower and upper boundary. 16 Almost the same line 
shapes are observed for 0.7 < A < 1. When we further decrease the value of A, a small peak 
near the upper boundary emerges for A < 0.5, and grows as A is further decreasing. In the 
limit A — > 0, this peak eventually becomes the singularity at the upper boundary, where the 
line shape coincides with the density of states. 
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Fig. 3. (a)Density of states for psinon-antipsinon excitations at q = ir/2 , M z = N/A from data for 
iV=4096. 

(b) Spectral- weight distribution of psinon-antipsinon excitations in S zz (q,cj) at q = n/2 , M z = 
N/4 from data for 7V=4096. 
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5. Conclusion 

We have evaluated the longitudinal dynamic spin structure factor S zz (q, to) for the massless 
(0 < A < 1) spin- 1/2 XXZ chain at q = it/ 2 with half the saturation magnetization M z = 
N/4, for the system size N = 4096. We have calculated the relative contribution of psinon- 
antipsinon excitations and confirmed that they are also dominant in the critical regime. We 
have also seen that the observed critical exponents of the singularity at the lower boundary 
agree with the prediction of conformal field theory. 

We have observed here, for the first time, that the A-dependece of transition rates and 
the spectral weight distribution in the magnetic field. It can be seen that the transition rates 
coverges toward the constant function with the value 1 in the limit A — ► 0. An interesting 
behavior of line shapes for S zz (it/2,uj) is observed around 0.1 < A < 0.5, where a small peak 
emerges near the upper boundary and grows as A — ► 0. We hope this peculiar behavior of the 
peak will be measured in inelastic neutron scattering experiments. 
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